Temperature dependence of QGP viscosity over entropy ratio from hydrodynamical 
analysis of ALICE data in yfs^=2.7Q TeV Pb+Pb collisions 
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Within Israel-Stewart's theory of dissipative hydrodynamics, we have analyzed ALICE data for the 
centrality dependence of charged particles multiplicity, elliptic flow and pr spectra in ^/sj^j^=2.76 
TeV Pb-fPb collisions and obtained the temperature dependence of the QGP viscosity over the 
entropy ratio {{^) ^gp)- If temperature dependence of (f is parameterized as (7)^^^ = a-^-jr^ + 
experimental data favor a in the range 0-0.2. a > 0.4 is not favored by the data. Experimental 
data in ^/sj^j^— 200 GeV Au+Au collisions however, prefer a— OA. 

PACS numbers: 47.75.+i, 25.75.-q, 25.75.Ld 



I. INTRODUCTION 



In recent years there is considerable interest in QGP 
viscosity. Lattice simulations of QCD predicted that 
the strongly interacting nuclear matter undergoes a 
confinement-deconfinement cross-over transition [H-Q 
above a critical temperature Tc ~ 170 MeV. ^5^^=200 
GeV Au-|-Au collisions at Relativistic Heavy Ion Gol- 
lider (RHIC) [SnS,] and ys]^=2.76 TeV Pb-hPb colli- 
sions at Large Hadron Collider (LHC) [9l-[T2j. produced 
convincing evidences that in RHIG and LHG collisions 
a deconfined medium or Quark-Gluon Plasma (QGP) 
is produced. It is then important to characterize QGP 
in terms of its transport coefficients, e.g. shear viscos- 
ity, bulk viscosity and conductivity. Present paradigm 
is that shear viscosity over entropy ratio of any matter 
has a minimum, possibl y w ith a cusp, around the criti- 
cal temperature T = Tc [1^ . String theory based models 
(AdS/GFT) give a lower bound on viscosity of any matter 
v/s ^ l/47r [l4l . [Tsj . QGP viscosity over entropy ratio 
{{^)qgp) is largely uncertain. In a perturbative QCD, 
Arnold et al ^ estimated (^^^ ^ 1. In a SU(3) gauge 
theory, Meyer [l7| gave the upper bound (f)^^^ < 1; 
and his best estimate is ry/s=0.134(33)at T = 1.165Tc. 
At RHIC region, Nakamura and Sakai [I3I estimated the 
viscosity of a hot gluon gas as ?7/s=0. 1-0.4. In Figll] 
Nakamura and Sakai's 1^1 estimate of (-) is shown. 

^ ^ \ s / qgp 

The predictions have large uncertainty and it is difficult 
to conclude about the temperature dependence of QGP 
viscosity. 

In hydrodynamic models, several authors have esti- 
mated viscosity over the entropy ratio of the fluid pro- 
duced in Au-l-Au collisions at RHIC and in Pb-|-Pb col- 
lisions at LHC 19-28]. However, most of the studies 
ignore the temperature dependence of viscosity over en- 
tropy ratio. Recently, in [23,[33|, temperature dependent 
viscosity over entropy ratio was considered. Parameter- 



ized lattice QCD results ^Ts*] for {■^) was used. It was 
found that the elliptic flow in -/s^^ — 200 GeV Au-|-Au 
collisions at RHIC is dominated by the viscosity in the 
hadronic phase. It is largely insensitive to QGP viscos- 
ity. QGP viscosity dominates only at LHC energy -v/satjv 
— 2.76 TeV. Temperature dependence of QGP viscosity 
was studied in [27|, [11] . It was concluded that the tem- 
perature dependence of QGP viscosity over entropy ratio 
cannot be constrained by fitting spectra and elliptic flow 
data alone. 

Compared to Au-|-Au collisions at RHIC, in Pb-|-Pb 
collisions at LHC, fluid is produced at higher tempera- 
ture. QGP phase is much longer at LHC than at RHIC 
energy collisions. Temperature dependence of QGP vis- 
cosity over entropy ratio is then better determined from 
ALICE data. In the present paper, we analyze the AL- 
ICE data on the centrality dependence of charged parti- 
cles multiplicity, integrated and differential elliptic flow 
[i-[l3 in Vsa,jv=2-76 TeV Pb-|-Pb colhsions and deter- 
mine the temperature dependence of QGP viscosity over 
entropy ratio. Fluid is produced progressively at lower 
temperature as the collisions become more and more pe- 
ripheral. Centrality dependence of charged particles mul- 
tiplicity or elliptic flow will have the imprint of initial 
temperature and consequently will be sensitive to the 
temperature dependence of viscosity. Our strategy is 
simple, we examine whether a given temperature depen- 
dent (-) is consistent with the experimental data. 
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II. HYDRODYNAMIC EQUATIONS, 
EQUATION OF STATE, INITIAL CONDITIONS 
AND TEMPERATURE DEPENDENCE OF f^) 



We assume that in -^5^^^=2.76 TeV, Pb-|-Pb collisions 
a baryon free fluid is formed. The space-time evolution 
of the fluid is obtained by solving Israel-Stewart's 2nd 
order hydrodynamic equations. 



2 



0, 



(1) 



(2) 



Eq[T] is the conservation equation 
momentum tensor, T'"^ = {e + p)u^u'^ - 



for the energy- 



£, p 



and M being the energy density, pressure and fluid veloc- 
ity respectively. Eq|2] is the relaxation equation for the 
shear stress tensor tt'^". In Eql2l D = u^df^ is the con- 
vective time derivative, V<''w''> = ^(V^u'' + Vu'') - 
— uf^u^) is a symmetric traceless tensor, rj 
is the shear viscosity and is the relaxation time. It 
may be mentioned that in a conformally symmetric fluid 
relaxation equation can contain additional terms [sij . 
Assuming boost-invariance, Eqs[T] and [5] are solved in 
(r = \/t'^ — z'^,x, y, rjs = ^Inf^f) coordinates, with the 
code '"AZHYDRO-KOLKATA"', developed at the Cy- 
clotron Centre, Kolkata. Details of the code can be found 



32| 



Hydrodynamic equations are closed with an equation 
of state (EoS) p = p{e). In the present study, we 
use an equation of state where the Wuppertal-Budapest 
[sl lattice simulations for the deconflned phase is smoothly 
joined at T = Tc = 174 MeV, with hadronic reso- 
nance gas EoS comprising all the resonances below mass 
mres='2-5 GeV. Details of the EoS can be found in ^2^. 

Solution of partial differential equations fEas lll2p re- 
quires energy momentum tensor T'^'^ and shear stress ten- 
sor tt'^'^ at the transverse plane at the initial time n . We 
assume that in impact parameter b collision, at the ini- 
tial time Ti=0.6 fm, energy density is proportional to the 
transverse profile of the average number of binary colli- 
sions Ncoii(h,x,y), calculated in a Glauber model. 



^(b, X, y) = EiNcoiiih, X, y)]. 



(3) 



Si is the central energy density in a zero impact pa- 
rameter collisions. We also assume zero initial fluid ve- 
locity, Vx{x,y) — Vy{x,y) = 0. In dissipative hydrody- 
namics, apart from the energy density and fluid veloc- 
ity, shear stress tensor has to be initialised. We initialise 
the shear stress tensors to the boost-invariant values, 
^xx _ ^yy _ 2r]/3Ti, 77^^=0, T] being the shear viscosity 
coefficient. We have to specify rj as well as the relaxation 
time Tjr. For the relaxation time we use the Boltzmann 
estimate — 3r]/2p. Shear viscosity coefficient is param- 
eterized as follows: we assume that QGP viscosity over 
entropy ratio (^)ggp has a minimum at T = Tc and to 
be consistent with AdS/CFT conjecture we take I/Att as 
the minimum value. We approximate the temperature 
dependence of viscosity over entropy ratio of the QGP 
fluid by a linear relation. 
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FIG. 1: Three forms of temperature dependence of viscosity 
over the entropy ratio ((f)^^^), used in the present simu- 
lations are shown. The black circles and squares are lattice 
QCD predictions [3 for the QGP viscosity over entropy ratio. 
The while circles are approximate temperature dependence of 
(2) used in [H. 



the parameter a control the temperature dependence. 
Recognizing the uncertainty in the temperature depen- 
dence, we have considered three values for a, (i) stiff 
temperature dependence: a—QA (approximately consis- 
tent with lattice simulations), (ii) moderate tempera- 
ture dependence: a=0.2 and (iii) temperature indepen- 
dence: a—0. Below the critical temperature Tc, QGP 
fluid crosses over to hadronic fluid. For the hadronic 
fluid viscosity over entropy ratio {^)hrg use the re- 
sults for hadronic resonance gas [SSJ . The two viscosities 
are smoothly joined at T ^ Tc, 
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Particle production is not very sensitive to the 
smoothening parameter AT. In the following we use 
AT=0.02 GeV. In Fig[Tl temperature dependence of ^ 
for the three values of a is shown. 

We use the Cooper-Frye prescription to obtain the in- 
variant distributions of hadrons, from the freeze-out sur- 
face at Tf=130 MeV. In Cooper-Frye, invariant distribu- 
tion is obtained as. 



(6) 



where is the freeze-out hypersurface and f(x,p) is the 
distribution function. In viscous hydrodynamics, f{x,p) 
is a non-equilibrium distribution function. Assuming 
that the fluid is not far from equilibrium, we approxi- 
mate J{x,p) as , 



f{x,p) = feq{x,p)[l + (t)[x,p)], 



(7) 
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FIG. 2: The black circles are ALICE data for centrality de- 
pendence of charged particles multiplicity per participant pair 
in x/sjv]v= 2-76 TeV Pb+Pb collisions [lO]. The solid, long 
dashed, short dashed lines are hydrodynamic simulations with 
the three forms oi {-\ (see Eq|4|. 



where feq{x,p) is the equilibrium distribution function. 
The deviation function (j){x,p) is approximated as, 

where, e, P and T are (local) energy density, pressure 
and temperature respectively. 

We may mention here that there is some overlap be- 
tween the main idea in the present paper and that of 
Ref . [29! . We emphasize the differences here. In [29| . 
effect of temperature dependent viscosity over entropy 
ratio, on elliptic flow and transverse momentum spec- 
tra was investigated. A fixed temperature dependent 
('s)(;gp' fro™^ lattice QCD calculations, was used in the 
simulations. In the present simulations, on the other 
hand, we try to determine the temperature dependence 
of (f)qgp- Three different temperature dependence of 

(^) is tested against the experimental data. As shown 
in FiglH temperature dependence of (s)^^^ in is 
steeper than the steepest QGP viscosity over entropy ra- 
tio used presently. There are some minor differences in 
details of the simulation also, e.g. in equations of state, 
initialisation of shear stress tensor, decoupling tempera- 
ture etc. 



III. RESULTS 

A. Pb+Pb collisions at LHC 

The central energy density Si of the fluid, at the initial 
time Tj;=0.6 fm, is the only free parameter in the model. 
We fix £i to reproduce experimental charged particle's 
muhiphcity ^ = 1601 ± 60 in 0-5% Pb-|-Pb colhsion 
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FIG. 3: Black circles are charged particle's spectra in 0- 
5% Pb-fPb collisions The solid, long dashed, and short 
dashed lines are hydrodynamic model simulations for charged 
particle's spectra with the three forms of (^) (see Eq|4]). 



TABLE I: Central energy density (e;) and temperature (Ti) 
required to reproduce experimental charged particle's multi- 
plicity in 0-5% Pb+Pb coUisions at ^^^^=2.76 TeV at LHC. 
We have also listed Si and Ti (the bracketed quantities) re- 
quired to reproduce the experimental charged particle's mul- 
tiplicity in 0-5% Au+Au collisions at -ySiv]v~200 GeV at 
RHIC. Viscosity over the entropy ratio of the central fluid 
is also noted. 





{GeV/fm') 


(MeV) 


(-) 

\s J qgp 


a = 


135 + 7 


541 + 7 


0.08 




(32 + 2) 


(380 + 6) 


(0.08) 


a = 0.2 


78 + 4 


473 + 6 


0.43 




(22 + 2.0) 


(345 + 6) 


(0.20) 


a = 0.4 


49 + 2 


422 + 4 


0.96 




(16 + 1) 


(319 + 5) 


(0.34) 



[lo| . For the three forms of QGP viscosity, a—OA, 0.2 and 
0, we have simulated 0-5% Pb+Pb collisions and com- 
puted negative pion multiplicity. Resonance production 
is included. Noting that pion's constitute w80% of the to- 
tal charged particles, tt" multiplicity is multiplied by the 
factor 2 X 1.20 to compare with experimental charged par- 
ticle multiplicity. In tableU central energy density and 
temperature required to reproduce experimental multi- 
plicity in 0-5% collisions are listed. Uncertainty in Si or 
Ti reflect the uncertainty in ALICE measurements. As 
expected, central energy density or temperature of the 
fluid is reduced with increasing value of a. Charged par- 
ticles multiplicity is a measure of the final state entropy. 
In viscous fluid, during evolution, entropy is generated. 
Entropy generation increases with viscosity and more vis- 
cous fluid require less initial energy density to reach the 
fixed final state entropy. In tableHl we have also listed 
the viscosity over the entropy ratio of the central fluid. 
One note that central fluid viscosity is very large of if 
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FIG. 4: ALICE measurement ^13\ for charged particle's el- 
liptic flow in 4-particle cumulant method are compared with 
the simulated flows for three forms of (7)^^^ (see Eq|4]). 



^^\gp ^^^^ temperature dependence. 

However, stiff temperature dependence of ^ 
not consistent with the experimental centrality depen- 
dence of charged particles multiplicity. In Figl2] simu- 
lated charged particles multiplicity per participant pairs 
{-^ — ^^ch) for the three values of a, are compared 
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with the ALICE measurements [10]. ALICE data for the 
centrality dependence of charged particles multiplicity is 
not reproduced for a~OA. In mid-central and peripheral 
collisions, simulation produces fewer particles per partic- 
ipant pairs than in experiment. For moderate temper- 
ature dependence, a=0.2, agreement with experimental 
data in mid-central collisions improves but data in pe- 
ripheral collisions are still under predicted. Of the three 
forms, data are best explained for a=0, when (f )ggp is 
temperature independent. ALICE data on the central- 
ity dependence of charged particles multiplicity do not 
demand temperature dependent QGP viscosity over en- 
tropy ratio. Rather, it is consistent with the AdS/CFT 
minimal value (^) = 

\ S / qgp ilT 

Temperature independent (f )ggp a-lso best explains he 
ALICE measurements for the charged particle pT spectra 
in 0-5% Pb-fPb collisions [llj]. In Figl^l the simulated 
charged particles spectra, for the three values of a, a—Q, 
0.2 and 0.4, are plotted against the ALICE data. Note 
that all the three curves reproduce the charged parti- 
cles multiplicity (see FiglJ]). Spectra become flatter as a 
increases. ALICE data are best explained when (-) 

^ ^ s I qgp 

is temperature independent {a—Q). Even then, we do 
note that the model over predict the spectra at pT > 1 
GeV. The reason is understood. We have approximated 
charged particles spectra by the pion spectra. However, 
experimentally, kaon and proton spectra are flatter than 
the pion spectra. When constraint to reproduce the to- 
tal multiplicity, the approximation will overestimate the 
central fluid temperature, increasing the high pT yield. 
Integrated elliptic flow (U2) is an important observable 
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FIG. 5: The black circles are the ALICE measurement 
for charged particle's elliptic flow in 0-10%, 10-20%, 20-30% 
and 30-40% Pb-I-Pb collisions. The solid, long dashed and 
short dashed lines are hydrodynamical model simulations with 
the three forms of (f)^^^ (see Eq|4l). 



in relativistic energy collisions. It is a measure of col- 
lectivity in the medium. In FigHl ALICE measurements 
jl j | for V2 is shown. As the collisions become more and 
more peripheral, V2 or collectivity increases, till very pe- 
ripheral collisions where we observe a dip in collectivity. 
In FigHl solid, long dashed and short dashed lines are 
simulated (integrated) elliptic flow for q;=0, 0.2 and 0.4 
respectively. For q;=0.4, flow is largely under predicted. 
As it was for charged particles multiplicity, centrality de- 
pendence of elliptic flow also donor require stiff tempera- 
ture dependent {^) ^gp- Though the agreement with data 
is better for a=0.2, the model predictions under predict 
the data. Data are best described when a=0, i.e. (^) 
is temperature independent. 

A different conclusion is reached when we consider cen- 
trality dependence of differential elliptic flow {v2{pt))- In 
FiglSJ ALICE measurements for the charged particles dif- 
ferential elliptic flow in 10-20%, 20-30%, 30-40% and 40- 
50% Pb-f Pb collisions are compared with the simulated 
differential flows. ALICE data on differential flow also 
donor demand stiff temperature dependent (7)^^^- For 
a=0.4 simulated flows largely under predict experiment. 
Data are best explained when a^O.2, i.e. for moderate 
temperature dependence of (f)^^^- Temperature inde- 
pendent (f give comparatively poor description. 



B. Au+Au collisions at RHIC 



Present analysis of ALICE data in y/sj^j^=2.76 TeV 
Pb+Pb collisions clearly indicate that the data disfavor 
strong temperature dependent (7)^^^- A similar analy- 
sis can be performed for the experimental data in RHIC 
energy collisions. It is interesting to see whether or not 
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FIG. 6: (color online) In four panels, PHENIX data ^3^] for 
the charged particle's (differential) elliptic flow in 0-10%, 10- 
20%, 20-30% and 30-40% Au+Au collisions at ysjviv= 200 
GeV Au-|-Au collisions are shown. The solid, long dashed, 
short dashed lines in each panels are the hydrodynamic simu- 
lations for the pt spectra with the three forms of (2) (see 
Eqlll). 

QGP produced at RHIC also disfavors strong tempera- 
ture dependent (■^)ggp- As before, in •\/s^^=200 GeV 
Au-|-Au collisions also, we fix the central energy den- 
sity £i of the fluid to reproduce the charged particle's 
multiphcity ^ = 687 ± 37, in 0-5% Au-f Au collisions 
[3^ . They are listed in tablelH Compared to LHC energy, 
in RHIC energy collisions, fluid is produced at a lower 
temperature. For example, for temperature independent 
QGP viscosity (q!=0), the central fluid temperature is 
Ti «380 MeV, to be contrasted with T, «540 MeV in 
LHC energy collisions. For stiff temperature dependent 
QGP viscosity (a=0.4) the central fluid temperature is 
Ti Ki2>2Q MeV, approximately 100 MeV lower than the 
fluid temperature at LHC energy collisions. QGP pro- 
duced in RHIC energy collisions then explore more re- 
stricted temperature range than at LHC energy. Conse- 
quently, effect of viscosity will be comparatively less at 
RHIC than at LHC. 

In Figini PHENIX measurements for the charged 
particles (differential) elliptic flow in 0-10%, 10-20%, 20- 
30% and 30-40% Au-|-Au collisions are compared with 
the simulated flows. The solid, dashed and long dashed 
lines are present model simulations with the tempera- 
ture dependence controlling parameter a=0, 0.2 and 0.4 
respectively. As expected, elliptic flow is reduced as the 
temperature dependence of (f become stiffer. How- 
ever, reduction is much less than in Pb-|-Pb collisions at 
LHC energy. For example, in Pb-|-Pb collisions at LHC, 
in 10-20% collision, at pt=2 GeV, simulated flow is re- 
duced by ~ 20% when a is increased from to 0.2 and by 
~ 50% when a is increased from to 0.4. At RHIC energy 
however, the reduction is much less, ^6% and ~20% as a 
is increased from to 0.2 and to 0.4. Effect of viscosity 



is less on RHIC data due to reduced fluid temperature. 
Fits obtained to the experimental flow are also interest- 
ing. In central 0-10% collisions, for all the three forms of 
(■j)^^^, reproduces the experimental data upto px <1.5 

GeV. At higher pT data are under predicted. In 10-20% 
collisions, PHENIX data are marginally over predicted in 
simulation with temperature independent (7)^^^- Sim- 
ulations with moderate or stiff temperature dependent 
(^) better describe the data. In 20-30% or in 30-40% 

collisions also, temperature independent (7)^^^ over pre- 
dict the data at large px- Data are best explained with 
stiff temperature dependent (■j)^^^- We have not shown 
here, but our analysis of charged particles px spectra in 
Au-(-Au collisions at RHIC also favor stiff temperature 
dependent {^)ggpi rather than temperature independent 

or moderate temperature dependent {^)^gp ■ Unlike in 
Pb-I-Pb collisions at LHC, in Au-|-Au collisions at RHIC, 
experimental data prefer strong temperature dependent 

\s fqgp 

Present analysis of experimental data in RHIC and 
LHC energy collisions indicate that a unique temper- 
ature dependent (f)^^^ does not explains the experi- 
mental data both at RHIC and LHC energy. In LHC 
energy collisions, experimental data prefer temperature 
independent (a—0) or moderate temperature dependent 
(a=0.2) (fj^gpj but in RHIC energy collisions, data pre- 
fer strong temperature dependent (a=0.4) {^)ggp- Re- 
sults of the present analysis are similar to that of f27| . 
In [231, a hybrid model VISHNU was used to analyse 
the experimental data at RHIC and LHC energy colli- 
sions. In the hybrid model, below a switching temper- 
ature {Tswitch=^Q5 MeV), evolution of the fluid is gov- 
erned by the Israel-Stewart's 2nd order hydrodynamic 
equations. Above the switching temperature, the evolu- 
tion is obtained using the microscopic hadronic transport 
model from UrQMD. In UrQMD, the kinetic and chemi- 
cal freeze-out is governed by the Boltzmann equation and 
the model thus eliminates these two parameters required 
in pure hydrodynamic approach. In [27| it was concluded 
that the temperature dependence of (f cannot be de- 
termined from fitting spectra or elliptic fiow data alone. 
The present, pure hydrodynamic model analysis also in- 
dicate that an unique temperature dependent QGP vis- 
cosity over entropy ratio do not explains the experimental 
data both at RHIC and LHC energy. 



IV. SUMMARY AND CONCLUSIONS 

In Israel-Stewart's 2nd order hydrodynamics, we have 
analysed the ALICE data in ^j^j^=2.76 TeV Pb-^Pb 
collisions to determine the temperature dependence of 
QGP viscosity over entropy ratio. Parameterising the 
temperature dependence of QGP viscosity over entropy 
ratio as, (-) = a '^Z^" + for three values of a, 

' \s / qgp Ta iTT^ ' 
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a=0, 0.2 and 0.4, we have simulated Pb+Pb colhsions 
and compared with the ALICE data. ALICE data are 
best explained for a in the range 0-0.2. Strong temper- 
ature dependence a > 0.4 is disfavored by the data. In 
^/s^J^— 200 GeV Au-|-Au collisions, on the other hand, 
experimental data on elliptic flow or px spectra demand 
stiff temperature dependent (f)^^^- 

In the present analysis, we have assumed {^)^gp de- 
pends linearly on temperature. It is possible that the 
temperature dependence of {^)^gp is more complex. In- 
deed, one observes that for a=0-0.4, differential ellip- 
tic flow in Pb-I-Pb collisions at LHC, or in Au-f Au col- 
lisions at RHIC, are not properly explained in all the 
cohision centralities. While in ^5^^^=2.76 TeV Pb-|-Pb 
collisions a=0.4 largely underestimate V2{pt), a=0-0.2, 
overestimate V2{pt) at large pT- In y/sj^j^— 200 GeV 
Au-|-Au collisions also, elliptic flow at large px are over- 
estimated in peripheral collisions. Linear temperature 
dependence of (f may be inadequate for differential 
elliptic flow, which a very sensitive observable. Also, for 
hadronic gas viscosity over entropy ratio {'^) i^^^-, we have 
used the results of [SS^, it was not varied. Certainly, 
Ref.jSSj] does not provide the final word on how large 



r]/s is in the hadronic gas phase. Whether or not un- 
certainties in (-). . interferes with extraction of (-) 

V s / had \ s / qgp 

is not investigated. In hydrodynamical models, parti- 
cle production depends, apart from viscosity, also on a 
number of initial conditions, e.g. initial time, initial en- 
ergy density /velocity distribution, freeze-out condition 
etc. We have only varied the central energy density, did 
not explore all the possible initial conditions. Correla- 
tion between temperature dependence of (■^)qgp and var- 
ious initial conditions, e.g. CGC or fluctuating initial 
conditions etc. is not studied. Considering all the un- 
certainties, we summarize our analysis as follows: with 
Glauber model initialisation, if the temperature depen- 
dence of QGP viscosity over entropy ratio is parameter- 
ized as ( f ) = a + ^ , ALICE data in v^at ^ =2 . 76 
TeV Pb-|-Pb collisions are best explained for a in the 
range 0-0.2. Data disfavor strong temperature depen- 
dence, a > 0.4. However, RHIC data on charged par- 
ticle's elliptic flow or px spectra on the other hand pre- 
fer Q!=0.4, i.e. strong temperature dependent (f)^^^- 
A unique, linear temperature dependent QGP viscosity 
over entropy ratio fails to explain, simultaneously, the ex- 
perimental elliptic flow data at RHIC and LHC energy. 
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